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This talk is a rather sceptical review of our knowledge and understanding of deep 
inelastic scattering at low x (high energy). We show that the well known success of 
the DGLAP evolution equations in describing of experimental HERA data brought 
more problems and puzzles than answers. We advocate that more systematic 
theoretical and experimental investigations of the nonperturbative QCD are needed 
to clarify physics of DIS at high energy. 


1 Introduction 

In this talk we will answer two main questions: 

• What have we learned from deep inelastic scattering ( DIS ) at high 
energy ( low x ) from HERA? 

• What problems in DIS are still a challenge for QCD? 

Trying to find these answers we present here a critical review of our knowl¬ 
edge and understanding of DIS at high energy (low x). The motto of this talk 

n 

The well advocated success of the DGLAP evolution equation 0 in describ¬ 
ing HERA data in the region of low x brought more problems and puzzles than 
answers. To sort out these problems, we need to know more about nonpertur¬ 
bative QCD at small distances. 


2 Basics of QCD in DIS 

DIS occurs at small distances and this is the process most suitable to apply 
perturbative QCD (pQCD ). The advantage of DIS is the fact that we have 
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three general approaches: the renormalization group approach (RG), the Wil- 
sonJDperator Product Expansion ( WOPEEJ) and the factorization theorem ( 
FT0 ), these provide a deeper insight to the main properties of DIS, than any 
sophisticated calculation in pQCD. 

RG: Let us assume that we have integrated out all degrees of freedom 
with transverse momenta kt < A and obtained the effective Lagrangian: 

Left = Yl&WO*’ (!) 

i 

where A is the scale which separates the small and large momenta and Oi is 
arbitrary operator. The physical idea of RG is very simple, namely, physical 
observables ( dimensionless ) do not depend on scale A. This means that 
introducing a new scale A', we obtain a new effective Lagrangian with C l (A'). 
The RG says that C l ( A') = F{C l { A)) with known function F. Tins powerful 
method leads to the Dokshitser- Gribov - Lipatov - Altarelli - Parisild ( DGLAP 
) evolution equations that play the role of the Coulomb law in DIS. 

WOPE: This is a usual way to separate the power - like corrections (oc 
(^r)") from the logarithmic ones ( oc (aslnQ 2 )” ), where Q is the scale of 
hardness in our process ( the virtuality of the photon in DIS ). The WOPE 
maintains that any structure function (say F 2 ) can be written in the form: 

F 2 (x,Q 2 ) = F 2 LT (x,\nQ 2 ) + ± F 2 HT (x,\nQ 2 ) + O(^) , (2) 

where the cross section of the virtual photon is equal to a = ^ Q 2 m F 2 (x, Q 2 ). 
Both F 2 t and F^ T depend only on as In Q 2 and their In Q 2 dependence can 
be calculated in pQCD using the DGLAP evolution. The expansion of Eq.(2) 
is valid in any order of pQCD. 

FT: Allows us to separate the nonperturbative contribution (parton den¬ 
sities, f(x. A 2 ) ) at large distances (r t > from the perturbative ones ( 
coefficient functions, C(x, ln(Q 2 /A 2 )). For any structure function in Eq.(2) 
the FT gives: 

F 2 LT (x,Q 2 ) = Jdx'C(^M(Q 2 /A 2 )) ■ f(x',A 2 ) . (3) 

Coefficient function C can be calculated in pQCD while parton densities are 
the nonperturbative input in all our perturbative calculations. 

3 Success (!?) of the DGLAP evolution equations 
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3.1 What is the situation? 

The HERA experiment shows that the deep inelastic structure function F 2 (x, Q 2 ) 
has a steep behaviour in the small x region ( 10~ 2 > x > 10 “ 5 ), even for very 
small virtualities ( Q 2 ~ 1 GeV 2 ). Indeed, considering F 2 oc x~ x at Iowa;, the 
HERA data is consistent with a A which changes from 0.15 at C£ = 0.85 GeV 2 
to 0.4 at Q 2 = 20 GeV 2 . This steep behaviour is well describedO in framework 
of the DGLAP evolution equations, by all groups doing the global fit of the 
data ( GR\fl,MR90 and CTEQEj). No other ingredients such as the BFKL 
Pomeronla and/or the shadowing corrections (SCjllJ are needed to describe the 
experimental data starting from sufficiently low virtualities Q 2 « 1.5 GeV 2 . 
We now attempt to understand what compromise has been made to obtain a 
good description of the data and what has been actually done. 


3.2 What has been done? 


Let me recall a standard procedure of solving of the DGLAP evolution equa¬ 
tions. The first step: we introduce moments of the structure function, namely, 
xG(x , Q 2 ) = 2 F? /c e ~ U M(w, Q 2 ) duj , where contour C is located to the 

right of all singularities of moment M(lu,Q 2 ). The second step: we find the 
solution to the DGLAP equation for moment 


dM (lo,Q 2 ) 
d\nQ 2 


7(w)M(w,Q 2 ) . 


(4) 


The solution is 


M(uj,Q 2 ) = M(lo, Qq) ■ e 7G) MQ 2 /Qo) . (5) 

Here M(uj,Q g) is the nonperturbative input which should be taken from ex¬ 
perimental data or from “soft” phenomenology ( model). The third step: we 
find the solution for the parton structure function using the inverse transform, 
namely: 

xG(x,Q 2 ) = 

Therefore, to find a solution of the DGLAP equation we need to know the 
nonperturbative input M(w,( 5 g) and the anomalous dimension 7 (w), which 
we can calculate in perturbative QCD. We summarize below what has been 
calculated for the gluon anomalous dimension. We present the result as a table 
in which each element gives the order of the magnitude of the perturbative term 
that has been calculated. We use brackets (...) or [...] to indicate terms that 
have not yet been calculated. 


f ^ e u ln(l/x)+7(u>) ln(Q 2 /Qo) ^ 

Jc 2?r ® 


( 6 ) 
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Table 


We can now see what has been done in the global fits. The value of the 
anomalous dimension have been calculated in order ( two first rows in our 
table) and the nonperturbative input has been taken in the form M[us, Q q) oc 
UJ } UJo with Ldo ss 0.2 - 0.3. This means that the structure function at Q 2 = Q q 
increases as x _u10 at x —> 0 1. 

3.3 What is the price that must be paid? 

To understand the main feature of the low x behaviour of the deep inelastic 
structure functions, it is instructive to consider the asymptotics of xG(x,Q 2 ) 
using Eq.(6). This asymptotic is determined by the saddle point in lo = u>s 
and the equation for uss is 

ln(l/ x) + U=w s ln(Q 2 /Oo) = O . (7) 

Substituting 7 (w) = , one obtains us s = \J ~ - i.yi/^)^ anda;s 0 

in the region of low x. This means that us = wo contributes to the integral 
and, therefore, the energy behaviour of the DIS processes at high energy is 
determined by the initial partonic distributions-. For lower energies the saddle 
point contributes and the value of u>s ( see Ref. Ilil for the values of uss in HERA 
kinematic region) gives us the typical scale of us in the anomalous dimension as 
given in the Table. It turns out that in the HERAJdnematic region corrections 
of order of (^-) n should be important ( see Ref. Ell). 

3.4 The first lesson: 

We can describe the experimental data on the deep inelastic structure func¬ 
tion using only the DGLAP evolution equation without any new ingredients, 

“Strictly speaking this statement is correct for two global fits: MRS and CTEQ. The 
GRV fit has a different initial condition which we will discuss later. 
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provide we assume, that the initial parton distributions increases as and 
such behaviour has to be understood. The ordinary procedure of taking into 
account leading and next to leading contributions to the anomalous dimension 
cannot be justified since other corrections of the order of (^f)" with n > 2 
are large and should be included. An attempt to explain the steep initial dis¬ 
tribution was made in the GRV parameterization which starts from extremely 
low Q 2 ss 0.25 GeV 2 . It was shown that the DGLAP evolution can generate a 
power - like increase of the parton distribution at Q 2 « 1 GeV 2 . It is an inter¬ 
esting idea, but the open question is why, we can start the DGLAP evolution 
with so small virtualities where we have no reason to assume that only leading 
twist contribution is essential. We recall that the DGLAP evolution equations 
only describe the leading twist term in Eq.(2). 


4 A unified BFKL and DGLAP description of F\ data 


4-1 Low x resummation and the BFKL equation 

The contributions of the order of (^.) ra in 7 G (w) are given by the BFKL 
equations and were calculated in RefEd Namely, 


G ,N c OLS' ,N c OiSs r 

1bfkl( -) = > . M-) 


n= 1 


where 7 bfkl(^Zt) is given by iterations of 


( 8 ) 


1 = — x{1bfkl) where \(l) = 2^(1) - V’(t) - VK 1 - 1 ) • (9) 

71 LU 

The BFKL anomalous dimension of Eq.( 8 ) can be written in the form 






wr - to 

A 


( 10 ) 


where 7 bfkl^ 1 - 0 = wl) = \■ One can therefore see that the anomalous dimen¬ 
sion cannot exceed the value 7 = I or, in other word, the BFKL anomalous 
dimension should be essential in the kinematic region where 7 is close to 4. 


4-2 Where is the BFKL Pomeron? 

It is easy to find the kinematic region where the BFKL contribution become 
sizeable, since in the region of low x the solution of the evolution equation has 
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Contour plot for k and <y> GRV for Nucleon 



-1.0 0.0 1.0 2.0 3.0 4.0 

ln(Q 2 /GeV 2 ) 


Figure 1: Lines with fixed < 7 > and k for the GRV parametrization. 


a semiclassical form: 

xG(x,Q 2 ) = C • (i)<"> • (%) < ^ > , (11) 

x Qg 

where C,< to > and < 7 > are smooth functions of ln(l/x) and In (Q 2 /Qq). 
In Fig.l line < 7 >= \ is plotted, using the GRV parameterization, namely, 

< 7 > = dln ^(Q 2 / /g^ One can see, that HERA data has penetrated the 
region where the BFKL dymanics should be visible. On the other hand, as we 
have discussed, the data could be described without any contamination from 
BFKL. A possible answer has been given in two recent papers lilil, in which 
an attempt has been made to describe the HERA data using the following 
approach for 

7 g = + 1 bfkl [ Nc0 1) _ (12) 

7T U) 

In other word, the first row and all terms of the order ) n in the table has 

been included. Actually, such an approach was suggested long ago in Ref. 0, 
where the evolution equation which corresponds to Eq.(12) was written. The 
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main result of these two papers is as follows: using Eq.(12) we can describe 
all data starting with the flat (u>o —> 0) initial distribution at Q q « 1 GeV 2 . 
Such initial conditions give a natural matching with the “soft” processes. The 
recent breakthrough in the calculation of the next order correction to the BFKL 
equation E£l will lead to calculations in the nearest future of all terms in our 
table marked by [...]. It will allow us to calculate the next correction to Eq.(12) 
(Ay 6 ), namely 

AT 

A7 G = 7 “s + 7 osBfkl = + as y' _ ( 13 ) 

' 7TLU 
n— 1 

Using Eq.(13) we can evaluate the next order corrections and check how well 
our approach works. 

4-3 The second lesson: 

The widely held opinion that the DGLAP evolution equations describe the 
experimental data is too biased. The equations that incorporate the BFKL 
and DGLAP dynamics lead to better descriptions of the experimental data, 
providing an excellent matching between the “hard” and “soft” processes - flat 
initial parton distribution at Q 2 « lGeU 2 . 

5 Higher twists 

Everything, that has been discussed, so far is related only to the leading twist 
contribution to the deep inelastic structure function (first term in Eq.(2)). The 
only way, to separate the leading twist from the higher twist, is to consider 
sufficiently large value of Q 2 , where the higher twist contribution is expected to 
be small. However, practical estimates crucially depend on the Q 2 - dependence 
of F^ t . At first sight it appears that we do not know anything on F^ T . 
Actually, this is not true. We know a lot about the next order twist structure 
function: 

1. The physical meaning E3, namely,the high twist term is closely related 
to the correlation function of two partons in the parton cascade. 

2. The evolution equations!!-! These equations are similar to the Fadeev 
equations for many body problem, namely, they reduce the complicated parton 
interaction to the interaction of two partons with the same kernel as in the 
original DGLAP equations. 

3. The solution of the evolution equation in the region of very small x ( 
high energy) lla. 
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4. This solution suggests the following formula to fit the deep inelastic 
structure function: 


F 2 (x,Q 2 ) = F, 


DGLAP 


{x,Q 2 


m 

Q 


, 2 c(x)[F 2 DGLAP (x : Q 2 )f 


It should be stressed that this expression is quite different from that one 
which experimentalists assumed, namely, that F PT has the same Q 2 and x 
dependence as F 2 iT . 

However, nobody has tried to obtain numerical estimates for the high 
twist contribution, as the systematic computational approach to the evolution 
equations for F PT had not been-developed yet. 

This has recently been done £3 and the result is to some extent surprising. 
The extra power of Q 2 does not give the feeling that this contribution should 
be negligibelly small at least at Q 2 « lOGeH 2 . Indeed, if one wants to toi 

0 simple parameterization for F PT it is better to take F PT oc {F 2 T ) 2 ) £3 
and the Q 2 dependence due to anomalous dimension of F' 2 T compensates 
to a large extent for the 1 /Q 2 suppression. The importance of this result is 
obvious, since in all solutions of the DGLAP equations that there are on the 
market, it has been assumed that the higher twist contribution is small, and 
can be neglected at Q 2 = Qq , where Q q is the initial virtually of the photon 
from which we start the DGLAP evolution. Notice that in practice the value 
of Qq is rather small ( about 4 GeV 2 ). 


5.1 The third lesson: 

The common believe that the higher twist contributions are small at Q 2 ~ 
4GeH 2 does not look convincing. We have to study the higher twist con¬ 
tribution in detail as to develop a systematic computational approach to the 
evolution equations for the higher twist structure function. 


6 On the way to complete theory of DIS 

6.1 Lattice calculations 

For the first time the lattice calculations in DIS has achieved an accuracy 
that we have to discuss them seriously. We would like to recall, that the 
lattice experiment gives us solid information on the nonperturbative QCD 
contribution to DIS. In some sense this is the only way how we can obtain a 
selfconsistent and reliable nonperturbative contribution. 

The experimental errors in the lattice experiments are still large but, nev¬ 
ertheless, it gives a convincing result that the initial quark distributions at 




Figure 2: Reconstruction of u-quark distribution from lattice calculation ( dotted line) 
and the CTEQ parametrization (solid line) at fi 2 = 2.5 GeV 2 . 

Q 2 = 2.5 GeV 2 differs from experimental one ( see minireviewEl at DIS’97). 
Eig.2 shows the quark distribution predicted by lattice calculation ( see Ref. 
£il) and the same distribution that the CTEQ collaboration used as an initial 
one. Such a difference was expected since in the lattice experiment the leading 
twist contribution has been calculated while experimental data give the contri¬ 
bution of all twists at a definite value of the virtuality Q 2 . It is interesting to 
note that the leading twist quark distributions derived on the lattice turn out 
to be closer to one that was expected in the constituent quark model, where 
the mean momentum of the quark about i. 

6.2 The fourth lesson: 

The lesson from the first nonperturbative calculation is very simple: our usual 
method of separation leading twist from the higher twist does not work, the 
higher twist contribution is still large at Q 2 = Qq ss 3 — 4GeV 2 . 

6.3 A unique opportunity for a theoretical approach to DIS 

We would like to emphasize that the lattice calculations give an unique oppor¬ 
tunity to develop a self consistent theoretical approach to DIS. The following 
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strategy is advocated . First we use the lattice parton leading twist densities 
as the initial conditions for the GLAP evolution equations and solve them. 
The difference between experimental initial parton densities and the lattice 
one should be treated as the high twist contribution. For them we should use 
the high twist evolution equations,, and the theoretical status of the higher 
twist contribution was discussed. The above procedure will provide a theoret¬ 
ical approach to DIS and after it has been achieved we will be able to discuss 
DIS on the solid theoretical basis. However, the experts in lattice experiments 
still need to calculate the initial gluon density. 

7 Shadowing corrections: miracle or reality? 

The HERA data on deep inelastic structure function lead to puzzling result. 
On one hand, they can be successfully described in the framework of the 
DGLAP evolution equations without any new ingredients like the BFKL equa¬ 
tion and/or the shadowing corrections (SC),as have been discussed above. On 
the other hand, the parameter (k) which gives the estimate for the strength 
of the SC turns out to be large (n > 1) in theHSERA kinematic region (see 
Fig.l). This parameter k was estimated in RefsEIEI and it is equal to 

« = xG(x,Q 2 ) = xG{x,Q 2 ) . (14) 

To understand what is going on, we have to develop a theoretical approach 
in which we can treat the region of k ~ 1 in DIS. It should be stressed that 
previous attempts to develop a theory for the SC El only had a guaranteed 
theoretical accuracy for small k a s 1. Two such approaches were 
discussed recently: in the first oneEl pQGD was used at the edge of its validity 
( a s K < 1), while in the second ( see RefsEa ) the new approach was developed 
in the kinematic region of high parton density QCD. 

In Ref. El a new evolution equation was derived which describes that each 
parton in_the parton cascade interacts with the target in Glauber - Mueller 
approach El. The results are the following: (i) k is the correct parameter 
that determines the strength of the SC; (ii) the SC to the gluon structure 
function are big even in the HERA kinematic region, but nevertheless the 
value of the shadowed gluon structure function is still within the experimental 
errors or, another way of putting it, the difference between the shadowed and 
nonshadowed gluon structure functions does not exceed the difference between 
the gluon structure functions in the different parameterizations such as the 
MRS,GRV and CTEQ ones; (iii) the SC to F 2 (x,Q 2 ) in HERA kinematic 
region are so small that they can be neglected; (iv) the SC enter the game 
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before the BFKL equation and, therefore, the BFKL Pomeron cannot be seen 
in the deep inelastic structure function since it is hidden under substantial SC 
( it is interesting to note that this result is seen just from Fig.l where it is 
shown that the SC become important in the kinematic region where < 7 > is 
still smaller than i); and (v) in the region of low x the asymptotic behaviour 

of xG(x,Q 2 ) is xG{x,Q 2 ) —> ln(l/s) lnln(l/a;). This means that the 

gluon density is not saturatedEH unlike for the GLR equation. 


The new approach has been developed based on the ide. 
cal gluon field in the region of high parton density ( see Refscn 


iLthe semiclassi- 
pi andEZI L_The 

physical problem has been pointed out a long ago ( see Ref. E3and ReftZj for 
updated review): at high energy (low x ) and /or for DIS with heavy nucleus 
we are dealing with the system of partons so dense that conventional methods 
of pQCD does not work. However, the typical distances are still small for DIS, 
and this fact results in weak correlations between partons, due to small the 
coupling constant of QCD. 

a ie revolutionary idea, suggested in Ref. En and developed in Refs. O 
, is: in the Bjorken frame for DIS we can replace the complex QCD 
interaction between parton in such a system by the interaction of a parton 
i with energy fraction 27 with the classical field created by all partons with 
energy fraction x bigger than 27 . Indeed, in leading log(l/x) approximation 
of QCD all parton with x > 27 live for a much longer time than parton i, 
therefore, they create a gluon field which only depends on their xlensity-|Using 
this idea_and Wilson renormalization group approach, in Refsea andc3 ( see 
also Refea for elucidating remarks), the effective Lagrangian was obtained. 

It has been demonstrated that this-effective Lagrangian correctly repro¬ 
duces the DGLAP evolution equations Eil and even the BFKL Pomeron EZl in 
the limit of a sufficiently weak gluon field. However, the main problem in 
matching the two approaches: one which we discussed above in the beginning 
of this section and this one, is still an open problem. 


8 Crazy ideas 

Unfortunately, only two ideas on high energy behaviour of DIS we could call 
as crazy enough to be intemsting. 

The first one ( see Refl23 )is a generalization of the renormalization group 
approach to the longitudinal degrees of freedom. The arguments are based 
on the kt - factorization EZI and on similarity between InQ 2 and Ins. The 
answer is, roughly speaking, the BFKL amplitude, but with running QCD 
coupling which depends on energy. Certainly, this answer does not contradict 
unitarity and, perhaps, even the experimental data, but, of course, it is in 
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strong contradiction to the BFKL approach, that has been discussed in section 
4. None of the experts have seen a diagram which survives at high energy and 
in which the running coupling constant depends on s, and so it makes the 
whole approach rather suspicious. On the other hand, this approach is sure to 
stimulate a more detailed study of this problem. 

The second idea is more physical and based on the successful Buchmuller 
- Haidt parameterization of all data on phojxm - prpton interaction: photo¬ 
production and DIS E3. Accordingly to Refc3, the best parameterization of 
the HERA data is not the solution of the DGLAP evolution equations but a 
simple formula: 

p f ^ , (Q 2 + Ql) , xo 

F 2 (x 1 Q ) = mlog - log— , (15) 

G o x 

with m = 0.364; Xq = 0.04; and Q g = 0.55 GeV 2 . Eq.(14) gives a correct limit 
for F 2 at Q 2 —> 0 and describes the experimental data on the phntoproduction. 
Such a parameterization leads to the gluon structure functions xG(x, Q 2 ) = 
3 log which energy (x) dependence is just the same_as the solution to a new 
evolution equation which takes into account the SC E3. The physical idea is 
really crazy and sounds as follows: everything has happened at x of the order 
of 10 -2 , at such values of x the parton densities have reached their asymptotic 
values and in the region of x < 10 -2 we see only the asymptotic behaviour of 
the parton densities at low x. 

9 Resume 

I hope that I have convinced you that the low x behaviour of the deep inelastic 
structure function is still an open problem and the success of the DGLAP 
evolution equation caused more problems than answers. 

I believe, that I gave enough arguments to claim: 

1. That corrections to the DGLAP anomalous dimension are not small in 
HERA kinematic region at low x ; 

2. That the higher twist terms are not small at low x even for Q 2 « 
lOGeK 2 ; 

3. That the SC are not small in HERA kinematic region. 

Therefore, on one hand, we have to add the so called BFKL anomalous 
dimension to the DGLAP one in the region of small x at HERA, but, on 
the other hand, the BFKL evolution equation can not be detected due to the 
background of the strong SC. The explanation of why such different pictures 
can describe the experimental data, is closely related to the lack of direct 
measurement of the gluon structure function. Despites the beautiful data on 
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at HERA we only know the value of the gluon structure function within 
100 % errors. 

The only way to decide between our theoretical approaches is to predict the 
behaviour of the other processes at high energy, such as diffractive production, 
inclusive cross section and correlations in DIS both for nucleon and nuclear 
targets. Howevere, we have to admit that the accuracy of our theoretical 
calculations in all other processes, is not as good as for the total cross section 
of DIS. 

The good news is that the success in lattice calculations allows us to a 
large extent to control the nonperturbative contribution to DIS. Based on the 
lattice calculation, we can develop a selfconsistent theoretical approach ( see 
section 6.3) which enable us to clarify some of our problems in near future. 

A special chapter of DIS is the “hard” processes with a nuclear target, 
which I did not touch upon here. They certainly will help us to have a better 
understanding of both the higher twist contributions and the SC, but this is a 
subject for a separate talk. 

Acknowledgments 

I am very grateful to my colleques at HEP department in the Tel Aviv uni¬ 
versity and at DESY theory group for hot and provocative discussions on the 
subject. Special thanks to Asher Gotsman for his comments on the manuscript. 


References 


1 . 

2 . 

3. 


4. 


5. 

6 . 

7. 

8 . 

9. 


K. Wilson: Phys. Rev. 179 (1969) 1444. 

J.C. Collins, D.E. Soper and G. Sterman :Nucl.Phys. B308 (1988) 833. 
V.N. Gribov and L.N. Lipatov:5cw. J. Nucl. Phys. 15 (1972) 438; L.N. 
Lipatov: Yad. Fiz. 20 (1974) 181; G. Altarelli and G. Parisi:Auc/. Phys. 
B126 (1977) 298; Yu.L. DokshitsenSW Phys. JETP 46 (1977) 641. 

HI collaboration: S. Aid et a/., Nucl. Phys. B470 (1996) 3; 

ZEUS collaboration: M. Derrick et al.,Z.Phys. C69 (1996) 607, C72 


(1996) 399. 

R.G.Roberts -.Plenary talk at DIS’97 , Chicago, April 1997, RAL-TR-97- 


024, |hep-ph/9706269| . 


M. Gluck, E. Reya and A. Vogt: Z.Phys. C67(1995)433. 

A.D. Martin, R.G. Roberts and W.J. Stirling: Phys.Lett. B306(1993)145. 
CTEQ Collaboration, H.L.Lai et al.: Phys.Rev. D51(1995) 4763. 

E.A. Kuraev, L.N. Lipatov and V.S. Fadin: Sov. Phys. JETP 45(1977) 
199 ; Ya.Ya. Balitskii and L.V. Lipatov: Sov.J. Nucl. Phys. 28 (1978) 


13 





822; L.N. Lipatov: Sov. Phys. JETP 63 (1986) 904. 

10. L.V. Gribov, E.M. Levin and M.G. Ryskin: Phys. Rep. 100 (1983) 1. 

11. R.K. Ellis, Z. Kunst and E. M. Levin: Nucl. Phys. B420(1994) 517. 

12. T. Jaroszewicz:P/i 2 /s. Lett.B116(1982)291. 

13. R.S.Thorne: hep - ph /9701241 and references therein. 

14. J. Kwiecinski, A.D. Martin and A.M. Stasto: hep - ph /9703445. 

15. V.Fadin: Talk at DIS’97 , Chicago, April 1997. 

16. R.K.Ellis,W. Furmanski and R. Petronzio: Nucl.Phys. B207 (1982)1, 
B212 (1983) 29. 

17. A.P. Bukhvostov,G.V. Frolov, L.N. Lipatov and E.A. Kuraev: Nucl. 
Phys. B258 (1985) 601. 

18. J. Bartels: Phys. Lett. B298 (1993) 204, Z. Phys. C60 (1993) 471; 
E.M. Levin, M.G. Ryskin and A.G. Shuvaev: Nucl. Phys. B387 (1992) 


589. 

19. J. Bartels: Talk at DIS’97 , Chicago, April 1997. 

20. G. Schierholz: minireview, Talk at DIS’97 , Chicago, April 1997. 

21. T.Weigl and L. Mankiewicz: Phys. Lett. B339 (1996) 334. 

22. A.H. Mueller and J. Qiu: Nucl. Phys. B268 (1986) 427. 

23. A.L. Ayala, M.B. Gay Ducati and E.M. Levin: TAUP 2432 - 97, hep - 
ph/9706347; Nucl. Phys. B493(1997) 305. 

24. J. Jalilian-Marian, A. Kovner, L. McLerran and H. Weigert: hep - ph 
/9606337; H. Weigert: Talk at DIS’97 , Chicago, April 1997. 

25. A.H.Mueller: Nucl. Phys. B335 (1990) 115. 

26. L. McLerran and R. Venugopalan: Phys. Rev. D49 (1994) 2233,3352, 
D50 (1994) 2225, D53 (1996) 458. 

27. J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert: hep - 
ph/9701284; A. Kovner: Talk at DIS’97 , Chicago, April 1997. 

28. E. Laenen and E. Levin: Ann. Rev. Nucl. Part. 44 (1994) 199. 

29. Yu.V. Kovchegov: Phys. Rev. D54 (1996) 5463, hep - ph/9701229; 
Yu.V. Kovchegov, A.H. Mueller and S. Wallon: hep - ph/9704369. 

30. R.D. Ball and S. Forte jhep-ph/9703417| ; R.D. Ball: Talk at DIS’97 , 
Chicago, April 1997. 

31. S.Catani,M. Ciafaloni and F. Hautmann: Phys. Lett. B242 (1990) 
97, Nucl. Phys. B366 (1991) 135; J.C. Collins and R.K. Ellis: Nucl. 
Phys. B360 (1991) 3; E.M. Levin,M.G. Ryskin,Yu.M.Shabelsky and 
A.G. Shuvaev: Sov. J.Nucl.Phys. 53 (1991) 657. 

32. W. Buchmueller and D. Haidt: DESY 96-061, D.Haidt: Talk at DIS’97, 
Chicago, April 1997. 


14 



